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1.0  SUMMARY 


The  primary  objective  of  this  study  is  to  understand  the  fundamental  mechanisms  of 
interacting  coated  bodies  through  the  tribological  modeling  and  to  develop  optimized  coating 
systems  for  severe  contact  applications.  Ultimately  the  tribological  modeling  provides  guidelines 
for  the  design  and  materials  development  of  tribological  coating  systems  in  the  future. 

The  emphasis  of  this  investigation  has  been  made  on  the  development  of  numerical 
models  to  study  the  stress  distribution  and  elastic  plastic  deformation  generated  inside  the 
multilayered  functionally  gradient  coating  systems  developed  for  the  heavy  load  and  scratch 
resistant  applications.  Finite  element  analysis  (FEA)  models  were  developed  for  the  two 
dimensional  and  three  dimensional  contact  problems.  A  widely  used  commercial  software 
ABAQUS/Standard  with  ABAQUS/CAE  pre  and  post  processors  was  used  for  the  development 
of  FEA  models.  The  project  was  divided  into  two  parts:  (i)  stress  modeling  and  (ii) 
nanoindentation  and  deformation  modeling. 

The  first  part  of  the  study  was  focused  on  the  study  of  stress  distribution  in  layered  and 
gradient  tribological  coating  system.  The  critical  parameters  such  as  coating  thickness, 
composition,  applied  load,  material  properties,  and  friction  were  identified  and  investigated  for 
their  effects  on  stress  distribution  inside  the  layered  and  gradient  tribological  coating  system. 

Aside  from  the  mechanical  stress  modeling  of  the  coating  system,  experimental 
nanoindentation  technique  in  connection  with  FEA  model  simulation  was  also  used  to  measure 
and  calculate  the  material  properties  such  as  hardness,  modulus  of  elasticity,  and  yield  strength  of 
the  coating  materials.  Therefore,  the  second  part  of  the  study  was  focused  on  the  material 
properties  measurement  and  calculation  by  the  aid  of  nanoindentation  experiment  in  connection 
with  FEA  model  simulation  and  the  study  of  elastic  plastic  deformation  inside  the  multilayered 
functionally  gradient  Ti/TiC/DLC  coating  system  using  the  measured  material  properties. 

2.0  INTRODUCTION 

When  two  bodies  are  in  contact  with  relative  motion  while  supporting  load,  stress  and 
deformation  fields  are  generated  below  the  contact  surface.  The  service  life  of  machine 
components  is  major  concern  to  the  users  and  it  is  greatly  influenced  by  the  friction  and  wear 
characteristics  of  the  contacting  surfaces.  Tribological  behavior  of  the  mating  surfaces  can  be 
significantly  improved  by  depositing  thin  hard  coatings  on  them.  Thin  layer  coatings  such  as 
TiN,  TiC,  CrN,  and  TiCN  deposited  on  the  contacting  surfaces  have  been  developed  for  various 
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applications  to  reduce  the  friction  and  wear  of  engineering  components  that  experience  severe 
operating  conditions.  Recently,  diamond  like  carbon  (DLC)  coatings  were  recognized  as  wear 
protective  coatings  since  amorphous  DLC  materials  have  the  unique  combination  of  high 
hardness,  low  friction  and  low  wear.  Although  DLC  coatings  have  the  superior  friction  and  wear 
characteristics  to  the  conventional  thin  layer  coatings,  DLC  coatings  deposited  on  soft  steel 
substrates  have  some  disadvantages  such  as  residual  compressive  stress  and  poor  adhesion 
strength. 

Lyubimov  et  al.  [1]  found  that  interfaces  between  layers  are  the  prime  locations  for  the 
film  failure,  e.g.  by  an  interface  crack  growth,  because  of  the  considerable  mechanical  property 
differences  at  the  interfaces.  Voevodin  et  al.  [2]  showed  that  DLC  coatings  on  steel  substrates 
have  the  application  limitations  such  as  contact  pressure  (below  1  GPa)  and  coating  thickness 
(below  0.5  pm)  in  order  to  avoid  the  coating  failure,  i.e.  brittle  fracture  and  delamination,  at  the 
coating/substrate  interfaces  due  to  significant  material  property  differences.  The  limitations  of 
DLC  coatings  were  overcome  by  applying  multilayer  design  concepts  to  the  coating  system. 
Voevodin  et  al.  [3]  proposed  the  several  design  concepts  for  the  architecture  of  multilayer 
coatings  to  improve  the  limitations  of  the  DLC  coatings.  They  suggested  the  functionally 
gradient  (FG)  Ti/TiC/DLC  coatings  increase  the  load  support  and  adhesion  of  the  DLC  coatings 
while  maintaining  low  friction  at  the  contacting  surfaces,  hi  such  coatings,  the  composition  and 
hardness  vary  gradually  from  the  substrate  to  DLC,  avoiding  drastic  variations  in  material 
properties.  In  addition,  Voevodin  et  al.  [4]  developed  the  optimized  design  of  hard  coatings 
using  the  multilayer  and  FG  approach  experimentally.  They  found  that  the  design  architecture 
provided  the  gradual  increase  of  the  hardness  from  the  substrate  to  DLC  coating  and  substantially 
increased  the  adhesion  of  DLC  coating. 

Intelligent  design  of  multilayered  and  FG  coating  systems  requires  a  thorough 
understanding  of  layer  architecture,  chemistry  and  mechanical  properties.  Lyubimov  et  al.  [1] 
concluded  that  the  optimal  choice  of  thickness,  combination,  and  deposition  conditions  of  layers 
should  consider  their  stress  state  to  prevent  failure  of  multilayered  coating  systems.  Therefore,  it 
is  of  significant  importance  to  analyze  the  stress  distribution  of  the  multilayered  and  gradient 
coatings  and  substrate  through  the  development  of  finite  element  analysis  models  to  identify  the 
critical  parameters  that  constrain  the  deformation  and  subsurface  stresses  of  coating  systems 
caused  by  normal  contact  conditions  and  to  understand  the  failure  mechanism  of  the  coating 
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systems.  In  addition,  to  develop  the  numerical  models  accurately,  the  material  properties  such  as 
hardness,  modulus  of  elasticity,  and  yield  strength  need  to  be  measured  and  calculated  for  each 
thin  coating  material.  Due  to  the  thickness  of  the  thin  coatings,  conventional  tensile  test  cannot 
be  performed  to  generate  the  stress  strain  relationship  for  each  coating  material.  One  of  the  most 
widely  used  methods  to  measure  the  material  properties  for  thin  coating  materials  is  the 
nanoindentation  developed  for  characterizing  thin  coatings  and  layers  mechanically  because 
nanoindentation  of  thin  coatings  provides  a  relatively  simple  approach  to  quantify  the  mechanical 
properties  of  thin  coatings.  Due  to  the  thickness  of  the  thin  coatings,  only  shallow  indentation 
depth  is  allowed  to  extract  the  measurement  data  from  the  thin  coating  systems.  From  the 
nanoindentation  of  coating  systems,  the  load  vs.  indentation  depth  curves  are  obtained  and 
analyzed  to  determine  properties  such  as  modulus  of  elasticity  and  hardness.  The  measured  and 
calculated  material  properties  are  used  as  inputs  to  the  finite  element  analysis  (FEA)  models  to 
study  stress  distribution  and  corresponding  the  elastic-plastic  deformation  behavior  of  the  coating 
systems. 

This  present  study  focuses  on  the  study  of  stress  distribution  and  elastic  plastic 
deformation  on  the  multilayered  and  functionally  gradient  coating  systems  developed  at  the  Air 
Force  Research  Laboratory.  Therefore,  the  elastic  plastic  finite  element  analysis  model  has  been 
developed  for  the  multilayered  and  functionally  gradient  coating  systems  due  to  the  nonlinear 
nature  of  the  coating  system  using  a  commercial  software  ABAQU  S/Standard.  The  FEA  model 
has  been  used  to  understand  the  mechanisms  involved  in  tribological  contacts  for  a  wide  range  of 
material  properties,  coating  thickness,  composition,  applied  load,  and  friction.  Also,  the  FEA 
model  has  been  used  to  determine  the  stress  distribution  in  the  multilayered  FG  tribological 
coating  system  and  examine  the  elastic  plastic  deformation  of  the  coating  systems  during 
nanoindentation.  Nanoindentation  experiments  were  conducted  to  measure  the  material 
properties  such  as  hardness  and  modulus  of  elasticity  for  each  coating  material.  In  addition, 
nanoindentation  FEA  simulation  in  connection  with  the  nanoindentation  experiment  was 
performed  to  calculate  the  yield  strength  of  each  coating  material.  The  ultimate  goal  of  this  study 
is  to  develop  optimized  multilayered  FG  coating  systems  for  high  adhesion  and  scratch  resistant 
applications.  Further,  the  study  is  to  provide  the  design  guidelines  for  coating  design  and 
material  development  in  the  future. 
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3.0  A  STUDY  OF  STRESS  DISTRIBUTION  IN  LAYERED  AND  GRADIENT 
TRIBOLOGICAL  COATINGS 

3.1  Introduction 

Since  the  early  1950s,  most  of  the  early  works  were  focused  on  developing  the 
mathematical  formulations  and  analytical  solutions  for  the  plain  strain  contact  problems  such  as 
the  elastic  deformation  of  the  layered  media.  Meijers  [5]  proposed  mathematical  formulations 
for  the  contact  problem  of  a  rigid  cylinder  pressed  on  an  elastic  layered  media  using  asymptotic 
solutions.  Satisfactory  solution  for  all  values  of  c/b  (half  width  of  contact  area/layer  thickness) 
and  for  Poisson’  ratio  between  0.0  and  0.5  was  obtained  for  thin  and  thick  layers.  In  1970s, 
numerical  analysis  for  the  plain  strain  contact  problems  was  performed  by  many  researchers  to 
avoid  the  mathematical  complexity.  Gupta  and  Walowit  [6]  developed  a  mathematic  formulation 
for  the  normal  contact  of  layered  elastic  solids  using  Fourier  transform  approach  of  Airy  stress 
function.  Numerical  solutions  of  contact  pressure  and  size  of  contact  zone  for  a  wide  range  of 
layer  thickness  were  determined.  Chiu  and  Hartnett  [7]  solved  three  dimensional  counterformal 
contact  problems  involving  layered  solids  using  the  generalized  Boussinesq  solution.  King  and 
O’Sullivan  [8]  developed  a  two  dimensional  quasi-static  model  to  determine  the  stress  field  in 
layered  elastic  half-space  under  combined  normal  and  sliding  cylindrical  contact.  O’Sullivan  and 
King  [9]  extended  the  two  dimensional  model  to  the  full  three  dimensional  model  to  determine 
the  subsurface  stresses  under  the  combined  normal  and  sliding  contact  on  a  elastic  half  space 
with  a  layer  for  different  material  properties  of  layers  and  substrate  system.  Cole  and  Sayles  [10] 
developed  a  two  dimensional  numerical  model  for  the  contact  of  layered  elastic  bodies  with  real 
rough  surfaces.  They  found  that  rough  body  layered  contact  has  considerable  influence  on 
contact  pressure  behavior  than  smooth  body  layered  contact.  Also,  Kuo  and  Keer  [11]  developed 
a  three  dimensional  model  to  consider  a  multilayered  half  space  under  combined  normal  and 
sliding  contact  by  a  spherical  indenter  using  Hankel  transforms.  They  discussed  the  stress 
components  for  various  coating  thicknesses.  Polonsky  and  Keer  [12]  applied  the  fast  Fourier 
transform  (FFT)  technique  with  the  multi-level  multi-summation  technique  to  the  elastic  layered 
contact  problem.  Chudoba  et  al.  [13]  developed  a  novel  methodology  based  on  the  combination 
of  theoretical  image  loads  method  and  indentation  experiments  for  the  evaluation  of  the  response 
of  coated  substrates  to  mechanical  contact.  They  showed  that  the  effects  of  the  addition  of 
intermediate  layer  to  the  coating/substrate  systems.  Although  many  researchers  developed  the 
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numerical  models  for  the  layered  medium  contact  problems,  their  models  were  confined  to  the 
elastic  contact  behavior  of  the  coating  system. 

As  an  alternative  numerical  approach,  finite  element  analysis  has  been  carried  by  a  few 
researchers  for  the  contact  problem  of  layered  surfaces.  Komvopoulos  [14]  developed  a  two 
dimensional  finite  element  model  to  investigate  the  thin  layer  effects  on  elastic  semi-infinite 
solid  compressed  by  a  rigid  body  and  found  the  general  solutions  for  the  subsurface  stress  and 
deformation  fields.  Also,  the  parametric  analysis  of  layer  thickness  relative  to  Hertzian  contact 
width  size,  friction  coefficient  and  stiffness  of  layer  were  performed.  Komvopoulos  [15] 
extended  the  elastic  layered  media  finite  element  model  to  the  elastic-plastic  contact  problem  of 
TiN/Ti  coating  system.  The  critical  parameters  such  as  layer  thickness,  material  properties  of 
layer  and  substrate  materials  and  indentation  depth  on  the  threshold  of  plasticity  of  the  elastic- 
plastic  layered  media  indented  by  a  rigid  cylinder  were  examined.  He  found  that  the  coating 
thickness  plays  an  important  role  on  the  initiation  and  propagation  of  plastic  zone  and  developed 
a  deformation  map  for  TiN/Ti  coating  systems.  Bhattacharya  and  Nix  [16]  developed  an  elastic 
and  plastic  model  to  study  the  deformation  associated  with  indentation  of  the  thin  film/substrate 
system  using  the  finite  element  method  and  found  the  effects  of  the  elastic  and  plastic  properties 
of  the  materials  on  the  hardness  of  the  film/substrate  system.  Also  the  elastic  compliance  of  the 
film/substrate  system  was  determined  as  a  function  of  the  indentation  depth  from  the  load- 
displacement  curve.  Sun  et  al.  [17]  developed  an  elastic-plastic  axisymmetric  finite  element 
model  to  investigate  the  effects  of  important  parameters  such  as  thickness  of  coatings  and 
material  properties  of  substrates  on  the  various  TiN  coating/substrate  systems.  They  found  that 
substrate  strength  and  layer  thickness  play  significant  roles  on  the  load  capacity  of  the  systems, 
development  of  plastic  zone,  and  subsequent  coating  failure.  Krai  and  Komvopoulos  [18] 
developed  a  three  dimensional  finite  element  model  to  determine  the  subsurface  stress  and  strain 
fields  due  to  sliding  contact  on  an  elastic-plastic  layered  medium.  They  found  the  effects  of  layer 
material  properties,  friction  coefficient,  and  applied  load  on  the  stresses  in  the  layer  and  the 
substrate.  Ye  and  Komvopoulos  [19]  developed  the  three  dimensional  finite  element  model  to 
investigate  the  effect  of  residual  stress  in  the  layered  media  due  to  normal  and  sliding  loading  and 
unloading  conditions.  They  showed  the  optimum  residual  stress  in  the  layer  is  a  function  of 
loading  type  and  friction  coefficient. 
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In  the  present  study,  a  three  dimensional  finite  element  model  has  been  developed  to 
simulate  the  elastic  plastic  indentation  process  of  the  layered  and  gradient  Ti/TiC  coating  system 
in  order  to  determine  the  stress  distribution  and  elastic-plastic  deformation  inside  the 
multilayered  and  functionally  gradient  coating  system.  Also,  a  parametric  study  was  performed 
to  understand  the  effect  of  critical  parameters  on  the  stress  distribution. 

3.2  Description  of  the  Layered  and  Gradient  Ti/TiC  Coating  System 

Interfaces  between  layers  are  prime  locations  for  coating  failure  since  considerably 
different  mechanical  properties  between  the  layers  exist.  Therefore,  multilayer  coatings  with 
gradually  changing  composition  (functionally  gradient)  are  highly  desirable  to  avoid  the  drastic 
changes  of  mechanical  properties.  Figure  3.1  shows  the  typical  layered  and  gradient  Ti/TiC/DLC 
coating  displaying  the  compositional  variation  of  the  Ti/TiC/DLC  coating  deposited  on  top  of  the 
440C  steel  substrate  [3],  Figure  3.2  depicts  the  indentation  on  the  multilayered  and  gradient 
Ti/TiC  coatings  deposited  on  the  semi-infinite  medium  with  a  rigid  spherical  indenter  used  for 
this  study.  Each  coating  surface  is  assumed  to  be  smooth.  The  multilayered  and  gradient  Ti/TiC 
coatings  are  assumed  to  be  a  series  of  perfectly  bonded  layers  with  unique  material  properties  and 
layer  thickness  although  the  gradient  coatings  don’t  have  distinct  layers  in  reality,  that  is,  they 
normally  have  continuous  gradual  change  of  Ti  and  C  composition.  The  a-Ti  layer  is  deposited 
on  the  440C  stainless  steel  substrate  as  a  bond  layer  in  order  to  enhance  the  adhesion  between  the 
substrate  and  TiC  coatings.  Then  coatings  from  Ti.90C.10  to  Ti.30C.70  with  variational  Ti/C 
composition  are  continuously  deposited  on  top  of  the  bond  layer.  Figure  3.3  illustrates  the 
schematic  of  functionally  gradient  (FG)  Ti/TiC  coating  design  used  in  this  study.  On  top  of  the 
440C  stainless  steel  substrate, 


Figure  3.1.  Multilayered  and  functionally  gradient  Ti/TiC/DLC  coating  system  that  shows 
the  gradual  composition  variation  across  the  thicknesses  of  the  layers. 
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Figure  3.2.  Schematic  of  indentation  with  a  rigid  spherical  indenter  on  multilayered  semi¬ 
infinite  medium. 


Material 


Hardness 


Elastic 

Modulus 


Thickness 


Poisson’s 

Ratio 


Ti.30C.70 

29  GPa 

370  GPa 

100  -  1000  nm 

0.20 

Ti.5oC.5o 

20  GPa 

290  GPa 

100  nm 

0.20 

Figure  3.3.  Schematic  of  functionally  gradient  Ti/TiC  coating  design  with  compositional 

and  mechanical  properties  variation. 

a-Ti  is  added  as  a  bond  layer  with  50nm.  On  top  of  the  a-Ti  layer,  Ti  90C.10  is  deposited  with  50 
nm  thickness.  As  the  TiC  composition  changes,  the  gradient  coating  is  generated.  Each  layer 
has  its  unique  material  properties  such  as  hardness,  elastic  modulus,  and  Poisson’s  ratio. 

Specific  thickness  for  each  layer  is  used  for  this  application.  Figure  3.3  consists  of  the  material 
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property  values  used  for  the  present  study.  The  detail  for  TiC  gradient  coating  procedure  is  well 
described  in  Voevodin  et  al.  [3].  Table  3.1  contains  the  contact  conditions  such  as  normal  loads, 
coefficients  of  friction,  thicknesses  of  Ti.30C.70  layer,  indenter  radius,  and  bond  layers  used  in  this 
analysis. 


Table  3.1.  Contact  conditions. 


Applied  Normal  Load  (W) 

4.4x1 0’4 -0.15  (N) 

Coefficient  of  Friction  (p) 

0.0  -  0.6 

Thickness  of  Ti.3oC.7o  (t) 

100-1000  (nm) 

Radius  of  Indenter  (R) 

0.0002  (m) 

Bond  Layers 

a-Ti,  Ti  Alloy,  Chromium 

3.3  ABAQUS  Finite  Element  Model 

The  contact  problem  described  in  Figure  3.2  can  be  modeled  with  the  finite  element  mesh 
shown  in  Figure  3.4.  The  finite  element  model  is  to  analyze  contact  stress  distribution  in  layered 
and  gradient  Ti/TiC  coating  system  generated  by  a  rigid  spherical  indenter.  The  commercial 
finite  element  software  package  ABAQUS™/Standard  was  used  to  model  the  coating  system 
under  the  plane-strain  assumption.  Due  to  the  symmetry  of  the  system,  only  a  quarter  of  the 
computational  domain  was  modeled.  To  reduce  the  computing  time,  the  outer  portion  of 
computational  domain  was  modeled  as  infinite  body  using  infinite  elements.  In  addition,  very 
fine  mesh  was  used  for  the  gradient  coatings  adjacent  to  the  coatings/rigid  indenter  interface  near 
the  contact  zone  in  order  to  determine  the  stress  distribution  and  deformation  field  accurately. 
Relatively  coarse  mesh  was  used  for  the  region  far  away  from  the  contact  zone  to  reduce  the 
computing  time.  Nearly  30,000  nodes  and  elements  were  used  to  model  the  entire  system.  Two 
types  of  elements,  eight-node  linear  brick  reduced  integration  element  (C3D8R)  and  eight-node 
linear  one-way  infinite  element  (CIN3D8),  were  used  for  the  coatings  and  semi-infinite  medium. 
The  coatings  were  modeled  by  assigning  different  material  properties  for  the  elements.  The  rigid 
spherical  indenter  is  represented  by  a  rigid  spherical  body  with  a  radius  of  0.2  pm. 

Coulomb  friction  was  used  to  model  the  contact  between  the  top  coating  surface  and  the 
rigid  indenter.  Due  to  the  perfect  bonding  assumptions,  the  continuous  displacement  at  the 
coating  interfaces  and  coating/substrate  interface  is  applied. 
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Figure  3.4.  Schematic  of  three  dimensional  finite  element  analysis  model  developed  for  the  nanoindentation  with  a  rigid 
spherical  indenter  tip  on  multilayered  and  functionally  gradient  coating  system. 


Since  the  strain  hardening  is  not  considered,  elastic  perfectly  plastic  stress-strain 
relationships  were  used  for  the  gradient  coatings  and  the  substrate.  The  onset  of  plastic  yield 
begins  by  the  von  Mises  yield  criterion  for  the  isotropic  materials.  Figure  3.4  illustrates  the  side 
and  perspective  view  of  the  finite  element  model  developed  for  this  study.  Figure  3.4 
demonstrates  the  von  Mises  stress  inside  the  coatings  and  substrate  generated  by  a  rigid  spherical 
indenter  on  the  plane  of  symmetry,  x-y  plane  and  x-y  and  y-z  planes,  respectively.  The 
computation  was  performed  using  a  workstation  that  consists  of  9  dual-processor  900  MHZ 
Itanium2’s  with  2  Gigs  of  RAM  each.  Please  note  that  the  residual  compressive  stresses 
associated  with  the  process  of  deposition  are  negligible  because  the  layered  and  gradient  coating 
thickness  is  a  few  micros. 

3.4  Results  and  Discussion 

In  order  to  validate  the  three  dimensional  finite  element  model  developed  for  this  study, 
the  elastic  Hertzian  point  contact  between  the  440C  stainless  steel  substrate  and  the  rigid 
spherical  indenter  was  performed.  Figure  3.5  (a)  shows  the  normalized  Hertzian  point  contact 
pressure  distribution  generated  by  a  rigid  spherical  indenter  with  the  coefficient  of  friction, 
p  =  0.0.  Please  note  that  the  full  contact  pressure  distribution  was  generated  using  the  symmetry 
of  the  problem  although  a  quarter  of  the  computational  domain  was  modeled.  The  magnitude  of 
maximum  Hertzian  pressure  is  1.0  GPa  under  a  normal  load,  W  =  0.0035  N.  Figure  3.5  (b) 
depicts  the  contact  pressure  distributions  at  the  center  of  the  contact  along  the  axis  of  symmetry 
obtained  by  finite  element  and  analytical  solutions  for  maximum  Hertzian  contact  pressure,  1.0 
GPa  and  1.5  GPa,  respectively.  The  slight  discrepancy  between  the  finite  element  and  analytical 
results  at  the  contact  edge  is  due  to  the  node  discretization  size  of  the  contact  elements.  The 
results  indicate  that  there  is  good  agreement  between  the  finite  element  and  analytical  solutions. 

The  full  von  Mises  stress  distribution  on  the  plane  of  symmetry  was  plotted  using  the 
symmetry  of  the  problem.  Von  Mises  stress  is  taken  into  account  to  analyze  and  predict  the 
failure  in  the  contact  problems.  Figure  3.6  represents  the  von  Mises  stress  distributions 
generated  by  a  rigid  spherical  indenter  inside  the  440C  stainless  steel  substrate  on  the  vertical 
plane  of  symmetry,  i.e.,  X-Y  plane  for  a  normal  load 
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Dimensionless  Pressure  (P/R.) 


Figure  3.5.  Hertzian  contact  pressure  distribution  on  semi-infinite  medium  indented  by  a  rigid  spherical  indenter;  (a) 
normalized  contact  pressure  distribution  by  finite  element  solution  for  W  =  0.0035  N  and  (b)  contact  pressure 
distributions  at  the  center  of  contact  along  the  plane  of  symmetry  obtained  by  finite  element  and  analytical  solutions. 


X 


(a) 


Von  Mises  Stress/PH 


Figure  3.6.  Von  Mises  stress  distribution  in  semi-infinite  medium  indented  by  a  rigid 
spherical  indenter  for  normal  load,  W  =  0.0035  N  (p,  =  0.0);  (a)  perspective  view  of 
normalized  von  Mises  stress  on  the  plane  of  symmetry  and  (b)  contour  plot  of  normalized 

von  Mises  stress  on  the  plane  of  symmetry. 

W  =  0.003 5N  and  p  =0.0.  Figures  3.6(a)  depicts  the  perspective  view  of  normalized  von  Mises 
stress  at  the  center  of  the  contact  on  the  vertical  plane  of  symmetry.  This  figure  displays  the 
smooth  variation  of  von  Mises  stress  along  the  depth  direction,  Y,  as  well  as  along  the  surface 
direction,  X.  The  calculated  maximum  von  Mises  stress  is  0.635  Ph  located  approximately  at 
Y  =  0.5  below  the  surface.  Figure  3.6(b)  illustrates  the  contour  of  the  von  Mises  stress 
distribution  presented  in  Figure  3.6(a). 


12 


Figure  3.7  demonstrates  stress  component  profile  variations  below  the  surface  at  the 
center  of  the  contact  along  the  axis  of  depth  for  W  =  0.0035  N  and  p  =0.0.  Various  lines  are 
generated  by  the  finite  element  solutions.  Dots  are  generated  by  the  analytical  solutions  for  the 
same  contact  conditions.  ax  decreases  rapidly  as  the  depth  increases.  Due  to  the  symmetry,  az 
has  the  same  values  with  ax  along  the  depth  direction.  <jyhas  the  maximum  value,  1,  at  the 
surface  and  deceases  gradually  as  the  depth  increases.  amises  has  the  maximum  value  of  0.635  Ph 
located  approximately  at  Y  =  0.5  on  the  depth  axis.  In  general,  the  good  agreement  for  stress 
distribution  between  the  finite  element  and  analytical  solutions  was  achieved  at  the  various 
depths  below  the  contact  surface. 

Figure  3.8  shows  the  stress  component  profile  variations  below  the  surface  at  the  center 
of  the  contact  along  the  axis  of  depth  for  W  =0.0035  N  and  p  =0.0.  The  stress  components  such 
as  ax,  Gy,  gz,  and  amises  are  displayed.  The  diagram  depicts  the  layered  and  gradient  Ti/TiC 
coating  system  with  350  nm  Ti  30C.70  layer  thickness.  The  stress  component  profiles  display  the 
discontinuities  occurring  at  the  interfaces  due  to  the  material  property  differences. 
Fundamentally,  the  layers  behave  as  beams  bonded  to  the  substrate.  The  maximum  von  Mises 
stress  occurs  at  the  Ti.50C.50  and  Ti.70C.30  interface.  The  von  Mises  stress  difference  at  the 
interfaces  is  not  exceeding  0. 1  GPa.  gx  and  ctz  is  more  influenced  by  the  layered  and  gradient 
coatings  while  ay  is  less  influenced.  The  magnitude  and  location  of  maximum  von  Mises  stress 
can  be  altered  significantly  when  the  different  design  of  layered  and  gradient  coating  system  is 
used. 

Figure  3.9  represents  the  von  Mises  stress  distribution  in  the  layered  and  gradient  Ti/TiC 
coating  system  on  the  plane  of  symmetry  for  the  normal  loading  contact  conditions  with  350nm 
Ti.30C.70  layer  thickness,  W  =  0.0035  N  and  p  =0.0.  Figure  3.9(a)  shows  the  perspective  view  of 
the  von  Mises  stress  generated  by  a  rigid  spherical  indenter 
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Figure  3.7.  Analytical  and  finite  element  solutions  of  normalized  stresses  in  the  semi 
infinite  medium  due  to  the  normal  load  at  the  surface  (jj,  =  0). 
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Figure  3.8.  Analytical  and  finite  element  solutions  of  the  stresses  in  a  layered  and  gradient 
coating  system  due  to  normal  load  at  the  surface  (jo.  =  0.0). 
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(b) 

Figure  3.9.  Von  Mises  stress  distribution  in  layered  and  gradient  Ti/TiC  coating  system 
indented  by  a  rigid  spherical  indenter  for  indentation  load,  W  =  0.0035  N  (p,  =  0.0);  (a) 
perspective  view  of  normalized  von  Mises  stress  on  the  plane  of  symmetry  and  (b)  contour 
plot  of  normalized  von  Mises  stress  on  the  plane  of  symmetry. 

on  the  plane  of  symmetry.  Instead  of  having  smooth  von  Mises  stress  distribution  inside  the 
layered  and  gradient  Ti/TiC  coatings,  discontinuous  stress  distribution  at  the  layered  and  gradient 
coating  interfaces  exists  due  to  the  material  property  differences.  Overall  shape  of  von  Mises 
stress  distribution  is  similar  to  the  one  without  the  layered  and  gradient  Ti/TiC  coating.  The 
calculated  maximum  von  Mises  stress  is  795  MPa  located  at  Ti.soC  so  and  Ti.70C.30  interface 
along  the  depth  direction  below  the  surface.  Figure  3.9(b)  demonstrates  the  contour  of  the  von 
Mises  stress  distribution  presented  in  Figure  3.9(a). 

Figure  3.10  shows  the  effect  of  coating  design  on  von  Mises  stress  for  the  normal  load, 

W  =  0.0035  N  and  p=  0.0.  Figure  3.10  illustrates  the  von  Mises  stress  distribution  comparison 
between  a  700  nm  Ti.30C.70  layer  coating  on  440C  steel  substrate  and  a  700  nm  layered  and 
gradient  Ti/TiC  coating  on  440C  steel  substrate.  Please  note  that  the  von  Mises  stress 
distributions  are  obtained  at  the  center  of  the  contact  along  the  axis  of  depth  below  the  contact 
surface.  The  magnitude  of  von  Mises  stresses  occurred  inside  the  Ti.70C.30,  Ti.90C.10  and  a-Ti 
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coatings  are  smaller  than  the  one  of  the  monolayer  coating  because  of  the  low  moduli  of 
elasticity  for  Ti.70C.30,  Ti.90C.10  and  a-Ti  coatings.  Monolayer  coating  depicts  the  sharp  von 
Mises  stress  discontinuity  occurring  at  the  interface  due  to  the  high  modulus  difference  of  film 
and  substrate.  However,  the  layered  and  gradient  Ti/TiC  coating  represents  the  gradual  stepped 
von  Mises  stress  changes.  While  the  monolayer  coating  has  the  maximum  von  Mises  stress  at 
the  Ti.30C.70  and  substrate  interface,  the  layered  and  gradient  Ti/TiC  coating  has  its  maximum 
von  Mises  stress  at  the  Ti.50C.50/Ti.70C.30  interface.  The  magnitude  of  maximum  von  Mises  stress 
decreases  from  0.909  GPa  to  0.818  GPa  when  the  monolayer  coating  is  replaced  by  the  layered 
and  gradient  Ti/TiC  coating.  Therefore,  the  layered  and  gradient  coating  is  more  beneficial  than 
monolayer  coating  since  it  reduces  the  magnitude  of  maximum  von  Mises  stress  in  the  coating 
system.  In  addition,  it  results  in  lower  magnitude  differences  of  von  Mises  stress  at  the 
interfaces. 

In  order  to  represent  the  spatial  distribution  of  von  Mises  stress  in  the  layered  and 
gradient  Ti/TiC  coating  system,  the  results  from  the  FEA  model  are  presented.  Figure  3.1 1 
shows  contours  of  von  Mises  stress  in  the  layered  and  gradient  Ti/TiC  coating  system  indented 
by  a  rigid  spherical  indenter  on  the  plane  of  symmetry,  X  -  Y  plane. 


Figure  3.10.  Effect  of  coating  design  on  von  Mises  stress  in  (a)  monolayer  coating  system 
and  (b)  layered  and  gradient  coating  system  for  normal  load,  W  =  0.0035  N  (p,  =  0.0  ). 
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Figure  3.11.  Contour  plots  of  von  Mises  stress  in  layered  and  gradient  coating  system  indented  by  a  rigid  spherical  indenter 

for  normal  load;  (a)  W=  4.4  x  10  4  N,  (b)  0.0035  N,  (c)  0.0119  N,  and  (d)  0.0283  N. 


Figure  3.11  illustrates  that  the  von  Mises  stress  fields  below  the  surface  vary  as  the  load  increase 
from  W  =  4.4  x  10'4  N  through  0.0283  N  with  p  =  0.0.  Please  note  that  the  maximum  von  Mises 
stress  always  occurs  at  certain  depth  along  the  axis  of  depth  below  the  center  of  the  contact  due 
to  the  symmetry  of  the  problem.  Figure  3.1 1(a)  is  the  contour  plot  of  von  Mises  stress  for  a 
normal  load,  W  =  4.4  x  10"4  N.  The  maximum  von  Mises  stress  occurs  inside  the  Ti.30C.70 
coating  with  the  magnitude  of  415.7  MPa.  With  the  relatively  light  normal  load,  W  =  4.4  x  10‘4 
N,  only  shallow  indentation  depth  occurs  and  shows  the  elastic  behavior  of  the  coating/substrate 
system.  Figure  3.1 1(b)  is  the  contour  plot  of  von  Mises  stress  for  a  normal  load,  W  =  0.0035  N. 
The  maximum  von  Mises  stress  occurs  inside  the  Ti.50C.50  coating  with  the  magnitude  of  794.6 
MPa.  Figure  3.1 1(c)  is  the  contour  plot  of  von  Mises  stress  for  a  normal  load,  W  =  0.01 19  N. 
The  maximum  von  Mises  stress  occurs  inside  the  Ti.50C.50  coating  with  the  magnitude  of  1 .087 
GPa.  Figure  3. 1 1(d)  is  the  contour  plot  of  von  Mises  stress  for  a  normal  load,  W  =  0.0283  N. 
The  maximum  von  Mises  stress  occurs  inside  the  substrate  with  the  magnitude  of  1 .394  GPa.  In 
general,  as  the  applied  normal  load  increases,  the  location  of  maximum  von  Mises  stress  moves 
from  inside  the  layered  and  gradient  Ti/TiC  coatings  to  substrate.  Also,  as  load  increases,  wider 
half  Hertzian  contact  width  and  higher  von  Mises  stress  are  generated.  Due  to  the  material 
discontinuities  inside  the  layered  and  gradient  Ti/TiC  coatings,  there  is  discontinuous  stress 
distribution  at  the  interfaces  below  the  contact  surface. 

Figure  3.12  depicts  the  von  Mises  stress  distribution  on  the  plane  of  symmetry,  X-Y 
plane,  for  500nm  Ti.30C.70  layer  thickness,  W  =  4.4  x  10'4  N  and  p  =  0.0.  Figure  3.12(a)  is  the 
perspective  view  of  the  von  Mises  stress  distribution  generated  below  the  contact  surface  along 
the  axis  of  depth.  It  illustrates  that  the  maximum  von  Mises  stress  with  the  magnitude  of  0.408 
GPa  occurring  inside  the  top  Ti.30C.70  coating.  Due  to  the  discontinuous  material  properties  at 
the  coating  interface,  von  Mises  stress  distribution  displays  drastic  stress  variations  at  the 
interface  along  the  axis  of  depth.  Figure  3.12(b)  is  the  contour  plot  of  von  Mises  stress 
distribution  presented  in  Figure  3.12(a).  Figure  3.13  shows  the  effect  of  top  Ti.30C.70  coating 
thickness  on  von  Mises  stress  for  W  =  0.0035  N  and  p  =  0.0.  Depending  on  the  top  Ti.30C.70 
coating  thickness,  the  stress  distributions  along  the  axis  of  depth  at  the  center  of  the  contact  vary 
significantly.  With  1  OOnm  top  Ti.30C.70  coating  thickness,  the  maximum  von  Mises  stress,  628 
MPa,  occurs  at  the  Ti.90C.10  and  a  -  Ti  interface.  As  the  top  coating  thickness  increases  to 
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Figure  3.13.  Effect  of  top  Ti.30C.70  coating  thickness  on  von  Mises  stress  in  layered  and 
gradient  coating  system  for  various  Ti.30C.70  thicknesses  (W  =  0.0035  and  p  =  0.0). 

400nm,  the  location  of  maximum  von  Mises  stress,  818  MPa,  moves  to  the  Ti.50C.50  and  Ti.70C.30 
interface.  When  the  top  coating  thickness  is  700  nm,  the  maximum  von  Mises  stress,  875  MPa, 
occurs  at  the  Ti.30C.70  and  Ti.50C.50  interface.  Finally,  when  the  1000  nm  thickness  is  used,  the 
maximum  von  Mises  stress,  813  MPa,  occurs  inside  the  top  Ti.30C.70  coating.  The  thicker 
Ti.30C.70  coating  (1000  nm)  results  in  less  magnitude  of  maximum  von  Mises  stress  than 
intermediate  coatings  thicknesses  (400  nm  and  700  nm).  The  results  indicate  that  the  maximum 
von  Mises  stress  location  can  be  moved  vertically  along  the  axis  of  depth  depending  on  the 
coating  thickness  under  the  same  normal  load  while  the  von  Mises  stress  in  the  substrate  changes 
negligibly.  Also,  the  results  suggest  that  the  desirable  top  coating  thickness  changes  depending 
on  the  loading  conditions.  Therefore  the  top  layer  thickness  has  a  significant  effect  on  the 
magnitude  of  the  von  Mises  stress  developed  in  the  layered  and  gradient  Ti/TiC  coating  system. 

It  is  noted  that  the  maximum  thickness  of  the  top  Ti.30C.70  layer  is  affected  by  the  residual 
compressive  stresses  due  to  the  thermal  expansion  coefficient  differences. 

Figure  3.14  depicts  the  effect  of  top  coating  thickness  on  contact  pressure  distribution 
indented  by  a  rigid  spherical  indenter  on  the  contact  surface  under  W  =  0.0035  N  and  p  =  0.0. 
The  dashed  line  represents  the  pressure  distribution  generated  on  the  semi-infinite  medium 
without  coatings.  It  illustrates  the  maximum  Hertzian  pressure,  1 .0  GPa,  is  generated  under  the 
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normal  load,  W  =  0.0035  N.  When  the  layered  and  gradient  Ti/TiC  coating  system  is  applied,  the 
maximum  Hertzian  pressure  increases  due  to  the  stiffer  coating  layers.  The  dashpot  line  pressure 
distribution  is  generated  for  the  500nm  top  Ti.30C.70  layer  thickness  of  layered  and  gradient 
Ti/TiC  coating  system.  Due  to  the  stiffer  material  properties  of  the  coatings,  the  magnitude  of 
maximum  Hertzian  pressure  is  higher  with  the  narrower  half  contact  width  compared  with  the  no 
layer  Hertzian  contact  condition.  The  solid  line  pressure  distribution  is  generated  for  the  1  OOOnm 
top  Ti.30C.70  layer  thickness  of  layered  and  gradient  Ti/TiC  coating  system.  This  line 
demonstrates  highest  magnitude  of  von  Mises  stress  with  the  narrowest  half  contact  width  under 
the  normal  load,  W  =  0.0035  N.  In  general,  the  pressure  magnitude  increases  due  to  not  only  the 
increase  of  Young’s  modulus  but  also  the  increase  of  top  Ti.30C.70  coating  thickness  although  the 
same  normal  load,  W  =  0.0035  N,  is  applied.  In  addition,  the  half  Hertzian  contact  width 
decreases  as  the  top  coating  thickness  increases  because  the  same  magnitude  of  normal  load  is 
applied  to  the  layered  and  gradient  Ti/TiC  coating  system. 


Figure  3.14.  Effect  of  top  coating  thickness  on  contact  pressure  distribution  on  layered  and 
gradient  coating  system  indented  by  a  rigid  spherical  indenter  ( W  =  0.0035  N,  p  =  0.0 ). 
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Figure  3.15  shows  the  effects  of  applied  load  on  von  Mises  stress  in  layered  and  gradient 
Ti/TiC  coating  system  indented  by  a  rigid  spherical  indenter  for  various  normal  loads  with 
350nm  Ti.30C.70  layer  thickness  and  p  =  0.0.  The  magnitude  of  von  Mises  stress  in  the  layered 
and  gradient  Ti/TiC  coating  system  increases  as  the  applied  normal  load  increases  from  0.0004N 
to  0.0283  N.  The  maximum  von  Mises  stress  for  the  normal  load,  W  =  0.0004N,  is  416  MPa 
occurring  at  the  Ti.30C.70  and  Ti.50C.50  interface.  As  the  normal  load  increases,  the  magnitude  and 
location  of  the  maximum  von  Mises  stress  are  changed.  The  maximum  von  Mises  stress  for  the 
normal  load,  W  =  0.0035  N,  is  795  MPa  occurring  at  the  Ti.50C.50  and  Ti.70C.30  interface.  The 
maximum  von  Mises  stress  for  the  normal  load,  W  =  0.01 19  N,  is  1.09  GPa  at  the  Ti.50C.50  and 
Ti.70C.30  interface.  Finally  the  maximum  von  Mises  stress  for  the  normal  load,  W  =  0.0283  N,  is 
1 .39  GPa  at  the  Ti.50C.50  and  Ti.70C.30  interface.  Generally  the  magnitude  and  location  of 
maximum  von  Mises  stress  distribution  vary  significantly  depending  on  the  magnitude  of  normal 
load  when  the  layered  and  gradient  Ti/TiC  coating  system  is  used. 

Figure  3.16  depicts  the  effect  of  bond  layer  material  on  von  Mises  stress  for  W  =  0.0035 
N  and  p  =  0.0.  Bond  layer  materials  such  as  a-Ti,  Ti  alloy,  and  chromium  are  used.  Table  3.2 
contains  the  material  properties  such  as  modulus  of  elasticity  and  Poisson’s  ratio  for  three  bond 
layers.  There  are  significant  von  Mises  stress  changes  inside  the  bond  layer.  The  higher  the 
modulus  of  elasticity  is,  the  higher  the  von  Mises  stress  is.  The  negligible  von  Mises  stress 
changes  occur  in  the  other  coatings  and  substrate  even  if  the  substantial  von  Mises  stress 
variation  occurs  in  the  bond  layer.  Using  the  stiffer  bond  layer  like  chromium  increases  the 
resistance  to  initiation  of  plastic  deformation  in  the  bond  layer.  However,  it  creates  the  higher 
von  Mises  stress  differences  at  the  interfaces  above  and  below  the  bond  layer. 


Table  3.2.  Material  properties  for  bond  layers. 


Modulus  of  Elasticity  (GPa) 

Poisson’s  Ratio 

a-Ti 

140 

0.25 

Ti  Alloy 

103 

0.34 

Chromium 

279 

0.21 

In  order  to  investigate  the  friction  effect  at  the  contact  interface,  the  Coulomb  friction  was 
used  at  the  rigid  indenter  and  Ti.30C.70  coating  interface.  Figures  3.17  demonstrates  the  von 
Mises  stress  distribution  generated  in  350  nm  Ti.30C.70  layer  thickness  coating  system  at  the 
center  of  the  contact  along  the  axis  of  depth  under  the  normal  loading  conditions  for  W  =  4.4  x 
10‘4  N  and  W  =  0.0283  N.  Figure  3.17(a)  illustrates  that  the  effect  of  friction  coefficient 
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variation  from  p  =  0.0  to  p  =  0.3  on  von  Mises  stress  variation  near  the  surface  is  negligible 
under  the  low  loading  condition.  The  maximum  von  Mises  stresses  are  produced  at  the  Ti.30C.70 
and  Ti.50C.50  interface  for  p  =  0.0  through  0.3.  At  higher  loading,  W  =  0.0283  N,  there  are 
significant  von  Mises  stress  distribution  changes  near  the  contact  because  the  friction  (tangential) 
force  at  the  contact  is  proportional  to  applied  load  shown  in  Figure  3.17(b). 


Figure  3.15.  Effect  of  applied  load  on  von  Mises  stress  in  layered  and  gradient  coating 
system  indented  by  a  rigid  spherical  indenter  for  various  normal  loads 
(Ti.30C.70  =  350  nm,  p  =  0.0). 


Figure  3.16.  Effect  of  bond  layer  materials  on  von  Mises  stress  in  layered  and  gradient 
coating  system  for  W  =  0.0035  N  and  p  =  0.0. 


23 


Von  Mises  Stress  (Pa) 


Von  Mises  Stress  (Pa) 


Figure  3.17.  Effect  of  friction  coefficient  on  von  Mises  stress  in  layered  and  gradient 
coating  system  for  normal  load  (a)  W  =  4.4  x  10'4  N  and  (b)  W  =  0.0283  (Ti30C  70  =  350  nm). 


Also,  the  magnitude  of  maximum  von  Mises  stress  increases  as  the  friction  coefficients  increase. 
However,  the  location  of  von  Mises  stress  remains  the  same,  at  the  Ti.50C.50  and  Ti.70C.30 
interface. 

Figure  3.18  shows  the  effect  of  plasticity  on  von  Mises  stress  for  350  nm  Ti.30C.70  layer 
thickness,  W  =  0.15  N  and  (0.  =  0.0.  In  order  to  obtain  the  yield  strength  of  the  materials  used  for 
the  coating  system,  the  following  relationship  is  used  because  the  strain-hardening  is  negligible 
for  the  coating  materials  used  (Komvopolous,  [15]). 

_  hardness 

^  yield  ^ 

Linear  elastic-perfectly  plastic  stress-strain  law  was  applied  for  the  FEA  model.  Figure 
3.18(a)  demonstrates  that  the  contour  plot  of  von  Mises  stress  indented  by  a  rigid  spherical 
indenter  at  the  center  of  the  contact.  It  depicts  the  high  stress  concentration  inside  the  Ti.50C.50 
layer.  Also,  it  shows  that  the  maximum  von  Mises  stress  occurs  inside  the  substrate.  Figure 
3.1 8(b)  illustrates  von  Mises  stress  at  the  center  of  contact  along  the  axis  of  depth.  Significant 
plastic  deformation  is  developed  inside  the  Ti.30C.70  and  a-Ti  coatings  because  the  yield  strength 
of  the  materials  was  reached  as  the  normal  load  increases.  Therefore,  the  plastic  deformation  of 
the  coatings  is  initiated  at  the  bond  layer  interface.  Plastic  deformation  is  the  major  concern  in 
design  of  the  layered  and  gradient  coating/substrate  system  because  the  tribological  coating 
failure  usually  initiates  and  propagates  from  the  plastic  deformation  zone. 
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Figure  3.18.  Effect  of  plasticity  on  von  Mises  stress  in  layered  and  gradient  coating  system  for  Ti30C70=  350  nm, 
W  =  0.15  N,  and  p.  =  0.0;  (a)  contour  plot  of  von  Mises  stress  on  the  plane  of  symmetry  and  (b)  stress  profile  of  von 
Mises  stress  along  the  centerline  of  the  contact  for  elastic  material  properties  and  elastic-plastic  material  properties. 


3.5  Summary  and  Conclusion 

Stress  distributions  in  a  layered  and  gradient  TiC  coating/440C  steel  substrate  system 
indented  by  a  rigid  spherical  indenter  were  determined  using  the  finite  element  analysis  model. 
Design  superiority  of  the  layered  and  gradient  coating  to  monolayer  coating  was  proved  to 
improve  the  integrity  of  coating  system.  The  layered  and  gradient  Ti/TiC  coatings  generate 
smooth  discrete  stress  transition  at  the  coating  interfaces  when  the  chemistry  and  structure  of  the 
coating  is  varied.  Therefore,  the  stress  discontinuities  at  the  interfaces  are  important  design 
factor  to  consider.  Top  layer  thickness  is  very  influential  to  the  layered  and  gradient  Ti/TiC 
coating  system  because  it  changes  the  magnitude  and  location  of  maximum  von  Mises  stress 
distribution  inside  the  coating  system  significantly.  Also,  the  top  layer  thickness  changes  the 
magnitude  of  maximum  Hertzian  pressure  and  half  contact  width  under  the  same  load  condition 
due  to  the  combination  of  the  stiffness  and  thickness  of  the  layer  material.  In  addition, 
magnitude  of  the  applied  normal  load  affects  the  maximum  Hertzian  contact  pressure  and 
subsurface  von  Mises  stress  significantly.  Effect  of  friction  on  subsurface  von  Mises  stress  is 
small  for  low  loading  conditions.  However,  effect  of  friction  on  subsurface  von  Mises  stress  is 
substantial  for  high  loading  conditions.  Effect  of  the  bond  layer  materials  on  the  von  Mises 
stress  should  be  considered  because  yielding  in  the  layered  and  gradient  Ti/TiC  coating  system 
initiates  at  the  bond  layer  when  a-Ti  bond  layer  is  used. 

In  conclusion,  stress  distribution  in  the  layered  and  gradient  Ti/TiC  coating  system 
indented  by  a  rigid  spherical  indenter  depends  strongly  on  many  critical  parameters  such  as  the 
layer  thicknesses,  coating  design,  material  properties,  applied  load,  and  the  friction.  Especially, 
the  magnitude,  location,  and  discontinuity  of  the  von  Mises  stress  at  the  interfaces  play  important 
roles  on  the  interfacial  failure  of  the  coating  system.  Design  of  coating  systems  with  layered  and 
gradient  Ti/TiC  coatings  increases  the  limitation  of  their  applications  when  the  loading 
conditions  are  specified. 

4.0  NANOINDENTATION  AND  FINITE  ELEMENT  ANALYSIS  OF 

MULTILAYERED  AND  FUNCTIONALLY  GRADIENT  TRIBOLOGICAL 
COATINGS 

4.1  Introduction 

Nanoindentation  has  become  widely  used  for  the  characterization  of  coated  systems  over 
the  past  three  decades.  Doemer  and  Nix  [20]  developed  a  method  for  interpreting  hardness  and 
Young’s  modulus  from  the  data  obtained  from  depth-sensing  indentation  instruments.  Later, 
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Pharr  and  Oliver  [21]  proposed  the  improved  nanoindentation  technique  and  conducted  the 
experimental  measurements  of  thin  film  mechanical  properties  such  as  hardness  and  modulus  of 
elasticity,  that  is,  properties  are  derived  from  measurements  of  load,  displacement,  and  time. 

Page  and  Hainsworth  [22]  reviewed  and  summarized  the  abilities  and  limitations  of 
nanoindentation  techniques  to  characterize  the  mechanical  properties  of  thin  coatings.  They 
considered  these  nanoindention  systems  as  ultimate  mechanical  property  microprobes  that 
characterize  thin  coatings.  Nanoindentation  technique  has  been  used  by  Voevodin  et  al.  [2]  to 
evaluate  the  mechanical  properties  of  diamond-like  carbon  (DLC)  coatings  prepared  by  pulsed 
laser  deposition.  There  have  been  numerous  analytical  and  finite  element  models  developed  for 
the  nanoindentation  of  thin  coatings.  (Knapp  et  al.  [23,24,25],  Bouzakis  et  al.  [26],  Chudoba  et 
al.  [13]).  They  also  conducted  nanoindentation  experiements  for  various  types  of  coated  systems 
and  developed  the  procedures  for  obtaining  the  mechanical  properties  of  thin  coatings  based  on 
finite-element  modeling  of  nanoindentation.  They  emphasized  the  situation  when  substrate 
effects  for  thin  coating  mechanical  properties  measurement  are  dominant,  that  is,  the  thin  hard 
coatings  over  softer  substrates. 

In  this  present  study,  the  elastic  plastic  finite  element  analysis  model  using  measured 
material  properties  has  been  developed  and  used  to  determine  the  stress  distribution  in  the  FG 
Ti/TiC/DLC  multilayered  tribological  coating  system  and  examine  the  elastic  plastic  deformation 
of  the  coating  systems  during  indentation.  Due  to  nanoindentation  experiment  technique  in 
connection  with  nanoindentation  FEA  model  simulation,  yield  strengths  of  the  coating  materials 
were  calculated  using  the  measured  material  properties  such  as  hardness  and  modulus  of 
elasticity  for  each  coating  material. 

4.2  Functionally  Gradient  (FG)  Multilayered  Tribological  Coating 

Voevodin  et  al.  [3]  illustrated  the  need  for  the  functionally  gradient  (FG)  multilayered 
tribological  coatings  in  heavily  loaded  contact  applications  and  developed  the  multilayered  FG 
Ti/TiC/DLC  coating  system  shown  in  Figure  3.1.  Instead  of  having  diamond  like  carbon  (DLC) 
coating  on  top  of  the  substrate,  there  is  Ti/TiC  coating  system  that  has  gradual  Ti  and  C 
composition  changes  to  improve  the  load  support  and  adhesion  between  the  DLC  coating  and  the 
substrate. 
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4.2.1  Nanoindentation 

To  measure  the  material  properties  such  as  hardness  and  modulus  of  elasticity  for  FG 
Ti/TiC/DLC  coating  system,  a  MTS  nano  indenter®  XP  was  used  for  the  experiment  under  the 
carefully  controlled  environment.  Figure  3.2  depicts  the  schematic  of  nanoindentation  for 
multilayered  semi-infinite  medium  with  a  spherical  indenter  tip.  The  spherical  indenter  tip  is 
loaded  vertically  during  nanoindentation.  It  is  assumed  that  gradient  coatings  have  a  series  of 
perfectly  bonded  block  layers  with  unique  material  properties  and  layer  thickness.  However,  the 
gradient  coatings  don’t  have  distinct  layers  in  reality.  They  normally  have  continuous  gradual 
change  of  Ti  and  C  composition.  It  is  also  assumed  that  there  is  no  interfacial  stress  generated  by 
the  deposition  process  at  the  interfaces.  In  order  to  develop  an  accurate  nanoindentation  model 
for  the  FG  Ti/TiC/DLC  coating  system,  material  property  measurements  for  the  inputs  to  the 
FEA  model  are  required.  It  is  difficult  to  measure  the  material  properties  for  each  distinct  layer 
from  the  FG  multilayered  coating  system  because  the  material  property  measurements  for  thin 
coatings  require  very  sophisticated  instruments  due  to  the  magnitude  of  coating  thickness.  Nano 
indenters  with  nano-Newton  resolution  in  loading  and  nano-meter  resolution  in  displacement 
allow  very  accurate  measurements  of  the  load-displacement  curve  for  the  thin  coatings.  Figure 

4.1  demonstrates  a  portion  of  two  dimensional  axi-symmetric  FEA  model  of  the  nanoindentation 
simulation  using  a  conical  indenter  tip  for  coating  samples.  This  FEA  model  was  used  for 
various  coating  samples  to  calculate  their  yield  strengths  of  thin  coatings  using  measured 
hardness  (H)  and  modulus  of  elasticity  (E)  data.  Knapp  et  al.  [23]  developed  the  two 
dimensional  axi-symmetric  FEA  model  with  the  shape  of  the  conical  indenter  tip  to  verify  that 
the  conical  tip  with  a  certain  angle  has  the  same  area  versus  depth  relationship  for  the  Berkovich 
indenter  tip.  A  diamond  Berkovich  tip  was  used  for  the  nanoindentation  experiments  to  measure 
the  hardness  and  modulus  of  elasticity  of  thin  coatings.  Since  the  hardness  of  thin  coating 
materials  are  relatively  high,  the  deformation  of  the  diamond  indenter  tip  is  not  avoidable. 
Therefore,  the  deformable  indenter  tip  was  used  for  the  FEA  model.  Each  coating  sample 
consists  of  1  pm  thick  coating  and  440C  stainless  steel  substrate.  The  area  function  for  the 
Berkovich  indenter  tip  was  calibrated  through  a  series  of  nanoindentation  onto  fused  silica 
sample.  Nanoindentation  experiment  was  performed  for  the  carefully  prepared  samples  to 
minimize  the  measurement  errors. 
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Nano  Indenter  Tip 


Substrate  Thin  Film  Coatin9 

Figure  4.1.  Finite  element  analysis  (FEA)  model  of  the  nanoindentation  simulation  with  a 

conical  indenter  tip  for  the  coating  samples. 

Figure  4.2  illustrates  the  general  load  vs.  displacement  curves  generated  by  the 
nanoindentation  experiments.  The  red  dots  show  the  typical  load  vs.  displacement  curve  for 
elastic  deformation.  The  reversible  load  vs.  displacement  curve  was  generated  during  loading 
and  unloading  when  the  nanoindentation  causes  the  elastic  deformation  on  the  sample  surface. 
The  blue  dots  depict  the  typical  load  vs.  displacement  curve  for  elastic-plastic  deformation. 
During  loading,  the  indenter  tip  makes  a  permanent  impression  on  the  sample  surface  due  to  the 
plastic  response  of  the  sample.  The  load  vs.  displacement  curve  during  loading  provides  the 
hardness  information  of  the  materials.  During  unloading,  the  load  vs.  displacement  curve  is 
controlled  by  the  elastic  response  of  the  sample.  The  load  vs.  displacement  curve  during 
unloading  provides  the  modulus  of  elasticity  information  of  the  materials. 

From  the  nanoindentation  experiments,  the  hardness  and  modulus  of  elasticity  for  thin 
coating  materials  are  measured.  However,  to  develop  the  elastic-perfectly  plastic  FEA  model  for 
FG  Ti/TiC/DLC  coating  system,  the  yield  strengths  of  thin  coating  materials  should  be  obtained. 
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Although  there  are  no  measured  yield  strength  results  from  nanoindentation  experiments,  by  the 
aid  of  the  numerical  nanoindentation  simulation  using  a  FEA  model,  the  yield  strengths  of  thin 
coatings  are  able  to  be  calculated. 

•  •  •  •  Elastic  Deformation 

•  •  •  •  Elastic  Plastic  Deformation 


Figure  4.2.  Load  vs.  displacement  curves  for  elastic  deformation  and  elastic-plastic 
deformation  from  nanoindentation  experiment  respectively. 

Therefore,  interactive  procedure  between  FEA  model  simulations  and  nanoindentation 

experiments  is  required  to  calculate  the  yield  strengths  of  the  thin  coating  materials.  It  is 

assumed  that  the  plastic  deformation  of  a  coating  is  initiated  when  its  von  Mises  stress  generated 

inside  the  coating  becomes  greater  than  yield  strength  of  the  coating  material. 

4.3  Results  and  Discussion 

4.3.1  Yield  Strength  Calculation  Procedure 

The  procedure  to  calculate  yield  strength  of  a  coating  material  is  the  following.  The 

measurement  of  TiC/DLC  coating  material  properties  such  as  hardness  and  modulus  of  elasticity 

was  conducted  using  a  MTS  nano  indenter  with  a  Berkovich  indenter  tip.  The  measured 

modulus  of  elasticity  of  the  coating  material  is  used  as  an  input  to  the  FEA  model  for  the 

nanoindentation  simulation.  Also,  yield  strength  of  the  coating  material  should  be  used  as  an 

input  for  the  nanoindentation  simulation.  The  initial  yield  strength  of  the  coating  material  is 
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guessed  based  on  the  measured  hardness  value  of  the  material.  This  step  requires  an  iterative 
procedure  to  obtain  a  best-fit  load  vs.  displacement  curve  from  the  nanoindentation  simulation  by 
changing  the  yield  strength  input  value  for  the  simulation.  Therefore,  when  the  computational 
best-fit  load  vs.  displacement  curve  matches  the  experimental  averaged  load  vs.  displacement 
curve  closely,  the  guessed  yield  strength  value  is  considered  as  the  yield  strength  of  the  material. 
Yield  strength  of  each  sample  is  determined  from  the  nanoindentation  simulation. 

To  calibrate  the  indenter  tip  area  function,  a  series  of  indentation  was  performed  on  fused 
silica  (Si02)  because  SiC>2  is  a  standard  material  for  nanoindentation.  The  indentation  depth  was 
limited  to  100  nm  because  each  TiC/DLC  coating  sample  has  1  pm  coating  thickness.  To  avoid 
the  substrate  effect  on  the  material  property  measurement  of  thin  coating  materials,  1/10  rule  was 
applied.  Twenty  five  indentation  tests  were  performed  on  the  SiC>2  with  100  nm  indentation 
depth.  The  nanoindentation  experiment  results  for  SiC>2  are  in  shown  Figures  4.3  (a)  through  4.3 
(d).  Figure  4.3(a)  depicts  the  hardness  vs.  displacement  curves  for  SiC>2.  The  measured  average 
value  for  the  hardness  of  SiC>2  is  8.9  ±  0.4  GPa.  Figure  4.3(b)  demonstrates  the  modulus  of 
elasticity  vs.  displacement  curves  for  SiC>2.  The  measured  average  value  for  modulus  of 
elasticity  of  SiC>2  is  72.3  ±  3.5  GPa.  Hardness  values  were  taken  at  the  maximum  penetration 
depths,  and  moduli  of  elasticity  were  calculated  from  the  upper  unloading  portions  of  the  load- 
displacement  curves  (Voevodin  et  al.,  [3]).  Figure  4.3(c)  illustrates  the  load  vs.  displacement 
curves  of  Si02  corresponding  to  the  hardness  and  modulus  of  elasticity  curves  of  Si02.  During 
loading,  Si02  deforms  elastically  and  then  yields.  During  unloading,  the  elastic  deformation  of 
SiC>2  is  recovered  while  the  plastic  deformation  remains  as  a  permanent  impression  on  the  contact 
surface.  The  load  vs.  displacement  curve  comparison  between  the  experimental  results  and  FEA 
results  on  Si02  is  shown  in  Figure  4.3(d).  The  yield  strength  value  of  Si02  was  determined  by 
FEA  model  nanoindentation  simulation.  Since  the  measured  average  modulus  of  elasticity  of 
SiC>2  is  72.3  GPa  from  nanoindentation  test,  72.3  GPa  is  used  as  an  input  value  to  FEA  model 
simulation  for  Si02.  Circles  show  the  load  vs.  displacement  curve  using  average  nanoindentation 
experimental  results.  The  best-fit  solid  line  is  obtained  through  the  FEA  nanoindentation  model 
simulation.  Y  =  7.1  GPa  and  E  =  72.3  GPa  were  used  for  the  inputs  to  nanoindentation 
simulation  because  Y  =  7.1  ±0.3  GPa  and  E  =  73  ±  0.5  GPa  were  found  in  the  literature  (Knapp 
et  al.,  [25]).  The  results  indicate  that  the  best-fit  load  vs.  displacement  curve  from  FEA  model 
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Load  applied  on  Sample  (N)  Hardness  (GPa) 


Figure  4.3.  Nanoindentation  results  for  Si02;  (a)  hardness  vs.  displacement  plot,  (b)  modulus  of  elasticity  vs.  displacement 
plot,  (c)  load  vs.  displacement  plot,  and  (d)  load  vs.  displacement  plot  for  experimental  results  and  FEA  results. 


simulation  agrees  well  with  the  average  load  vs.  displacement  curve  from  the  nanoindentation 
experiment. 

Nanoindentation  experiments  were  conducted  for  440C  stainless  steel  because  440C 
stainless  steel  is  used  for  the  substrate  of  all  Ti/TiC/DLC  coating  samples.  Figure  4.4(a)  shows 
the  load  vs.  displacement  curves  for  440C  stainless  steel  during  loading  and  unloading  stages  of 
nanoindentation.  Due  to  the  surface  roughness  of  440C  stainless  steel,  experimental  curves  are 
scattered.  Because  of  the  nature  of  the  440C  stainless  steel,  higher  load,  ~1 .8  mN,  was  required 
to  make  indentations  upto  100  nm  in  depth  than  SiC>2.  Also,  the  slope  of  the  load  vs. 
displacement  curves  during  unloading  for  440C  stainless  steel  is  steeper  than  the  slope  of  the 
load  vs.  displacement  curves  during  unloading  for  SiC>2.  Figure  4.4(b)  displays  the  comparison 
of  load  vs.  displacement  curves  between  the  average  experimental  results  and  the  best-fit  FEA 
results.  The  average  experimental  load  vs.  displacement  curve  matches  the  best-fit 
computational  load  vs.  displacement  curve  when  Y  =  4.0  GPa  and  E  =  220  GPa  are  used. 
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Figure  4.4.  Nanoindentation  results  for  440C  stainless  steel;  (a)  load  vs.  displacement  plot  and  (b)  comparison  of 
load  vs.  displacement  plot  between  average  experimental  results  and  best-fit  simulation  results. 


4.3.2  Nanoindentation  for  Ti/TiC/DLC  Coating  Material  Properties 

TiC/DLC  coating  samples  were  prepared  by  a  hybrid  of  magnetron  sputtering  and  pulsed 
laser  deposition.  Crystalline  a-Ti,  TiC  and  amorphous  DLC  films  were  formed  at  100  °C 
substrate  temperature  by  varying  film  composition.  X-ray  photoelectron  spectroscopy  (XPS)  was 
performed  on  the  coatings  for  the  characterization  (Voevodin  et  al.,  [2]).  Figure  4.5  shows  the 
schematic  of  nanoindentation  on  coating  sample  to  measure  the  material  properties.  Vertical 
displacement  of  the  indenter  tip  is  applied  to  the  sample  during  the  nanoindentation.  The 
nanoindentation  process  was  simulated  by  moving  the  indenter  tip  vertically  with  controlled 
displacement  while  calculating  reaction  force  at  the  contacted  area.  440C  stainless  steel  is  used 
as  a  substrate  for  all  coating  samples  that  havel  pm  thick  coating  deposited.  For  the  simulation, 
the  coefficient  of  friction  between  the  indenter  tip  and  sample  surface  varies  from  0.0  to  1 .0. 
Table  4.1  consists  of  the  measured  material  properties  such  as  the  modulus  of  elasticity  and 
hardness  for  TiC/DLC  coating  samples.  The  results  indicate  that  the  hardness  and  modulus  of 
elasticity  changes  drastically  as  Ti  and  C  composition  changes. 

W 


1  pm 


Figure  4.5.  Schematic  of  nanoindentation  with  a  Berkovich  indenter  tip  on  coating  samples 
such  as  Ti,  Tii_xCx,  and  DLC  on  440C  stainless  steel  to  measure  the  material  properties. 

The  experimental  results  for  Ti.53C.47  coating  from  nanoindentation  are  shown  in  figures 
4.6(a)  through  4.6(d).  The  indentation  depth  to  Ti.53C.47  coating  sample  is  limited  to  100  nm  to 
avoid  the  substrate  effect  on  measured  material  properties.  The  figure  4.6(a)  shows  scattered 
data  for  hardness  vs.  displacement  curves  because  of  the  surface  roughness  effect  on  the 
measurement.  The  hardness  increases  as  the  nanoindentation  depth  increases  for  Ti.53C.47 
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coating.  The  average  hardness  of  Ti.53C.47  coating  from  the  nanoindentation  experiments  is  14.2 
±1.8  GPa.  The  thick  solid  line  represents  the  average  hardness  value  of  Ti.53C.47  coating  along 
the  indentation  depth.  The  figure  4.6(b)  depicts  the  modulus  of  elasticity  vs  displacement  curves 
for  Ti  53C.47  coating.  The  average  modulus  of  elasticity  of  Ti.53C.47  coating  is  231.1  ±  34.6  GPa. 
The  thick  solid  line  shows  the  average  value  of  the  modulus  of  elasticity  for  Ti.53C.47  coating 
along  the  indentation  depth.  The  figure  4.6(c)  illustrates  the  load  vs.  displacement  curves  for 
Table  4.1.  Mechanical  material  properties  such  as  modulus  of  elasticity  (E)  and  hardness 


(H)  for  TiC/DLC  coating  samples. 


Samples 

Material  Properties 

DLC  at  2x10''  Pa 

E  =  350  GPa,  H  =  33  GPa 

Ti.19C.8i 

E  =  230  GPa,  H  =  13  GPa 

Ti.28C.72 

E  =  280  GPa,  H  =  22  GPa 

Ti.53C.47 

E  =  240  GPa,  H  =  14  GPa 

Ti.90C.10 

E  =  145  GPa,  H  =  4.5  GPa 

440C  stainless  steel 

E  =  220  GPa,  H  =  9  GPa,  ay  =  4  GPa  (FEA) 

Ti.53C.47  coating  from  nanoindentation.  Due  to  the  roughness  of  surface,  the  data  is  scattered. 

The  figure  4.6(d)  shows  the  comparison  of  averaged  experimental  results  and  best-fit  FEA  results 
for  Ti  53C.47  coating.  The  average  load  vs.  displacement  curve  agrees  well  with  best-fit  load  vs. 
displacement  curve  when  Y  =  10  GPa  and  E  =  240  GPa  are  used  for  Ti.53C.47  coating 
nanoindentation  simulation. 
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Figure  4.6.  Nanoindentation  results  for  Ti  53C  47  coating;  (a)  hardness  vs.  displacement  plot,  (b)  modulus  of  elasticity  vs. 
displacement  plot,  (c)  load  vs.  displacement  plot  during  loading  and  unloading  stages,  and  (d)  comparison  of  measured  average 
load  vs.  displacement  curve  and  best-fit  load  vs.  displacement  curve  from  experiments  and  FEA  simulations,  respectively. 


Figure  4.7  demonstrates  the  schematic  of  functionally  gradient  (FG)  Ti/TiC/DLC  coating 
system  design  with  programmed  variation  of  compositional  and  mechanical  properties  proposed 
by  Voevodin  et  al.  [3].  Each  layer  has  its  unique  material  properties  and  thickness.  Hardness 
and  modulus  of  elasticity  of  Ti/TiC/DLC  coating  system  are  measured  from  coating  samples.  25 
indentation  locations  per  sample  were  averaged  to  obtain  the  statistical  values.  The  yield 
strength  for  each  coating  material  is  calculated  from  FEA  nanoindentation  simulation  using  the 
measured  hardness  and  modulus  of  elasticity.  These  are  the  values  used  for  inputs  to  FEA  model 
for  this  study.  The  material  properties  of  each  layer  were  linearly  interpolated  from  the  measured 
material  properties  data.  To  improve  adhesion  between  the  substrate  and  upper  carbide  layers,  a- 
Ti  is  used  as  a  bond  layer. 

Elastic  Yield  Pni*i^nn\ 

Material  Hardness  Modu|us  Strength  Thickness  P  ^ 


Figure  4.7.  Schematic  of  the  functionally  gradient  Ti/TiC/DLC  coating  design  with 
variation  of  composition,  dimension,  and  mechanical  properties. 

4.3.3  Finite  Element  Analysis  Model  Results 

Figures  4.8(a)  and  4.8(b)  show  a  schematic  of  nanoindentation  on  FG  multilayered 
coating  system  with  a  spherical  indenter  tip  and  its  corresponding  finite  element  analysis  model. 
Figure  4.8(b)  depicts  two  dimensional  axi-symmetric  finite  element  analysis  model  used  for  the 
nanoindentation  simulation  that  consists  of  loading  stage  and  unloading  stage.  The  radius  of  the 
diamond  spherical  indenter  tip  is  50  pm.  Since  the  hardness  of  the  DLC  coating  is  as  high  as  the 
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diamond  tip,  the  deformable  diamond  indenter  tip  is  modeled  for  the  nanoindentation  simulation. 
Although  the  nominal  dimension  of  the  diamond  indenter  is  50  pm  in  radius,  actual  dimension  of 
the  diamond  indenter  is  46.95  pm  in  radius.  ABAQUS/Standard  was  used  to  model  the 
nanoindentation  simulation  on  the  multilayered  FG  Ti/TiC/DLC  coating  system.  Due  to  the 
symmetry  of  the  system,  axi-symmetric  model  was  developed  to  save  the  computational  time. 
Also,  the  outer  portion  of  semi-infinite  medium  used  infinite  elements  to  reduce  the  computing 
time.  Two  different  types  of  elements,  four-node  bilinear,  reduced  integration  with  hourglass 
control  element  (CAX4R)  and  four-node  linear  one-way  infinite  element  (CINAX4),  were  used. 


(a)  (b) 

Figure  4.8.  (a)  Schematic  of  nanoindentation  on  multilayered  coating  system  with  a 
spherical  indenter  tip,  (b)two  dimensional  axi-symmetric  finite  element  analysis  model 
developed  for  nanoindentation  simulation  on  a  functionally  gradient  Ti/TiC/DLC 
multilayered  coating  system  with  a  spherical  indenter  tip. 

Figure  4.9(a)  shows  the  effect  of  multilayered  FG  Ti/TiC/DLC  coating  system  on  von 
Mises  Stress  distribution  during  the  nanoindentation  for  indenter  displacement,  6  =  50  nm  and 
friction  coefficient,  p  =  0.0.  Figure  4.9(b)  describes  the  schematic  of  FG  multilayered  coating 
system  design  used  for  this  study.  Figure  4.9(a)  illustrates  the  von  Mises  stress  profiles  along  the 
depth  direction  on  the  plane  of  symmetry  at  the  center  of  contact.  The  dotted  line  represents  the 
von  Mises  stress  profile  along  the  depth  direction  for  the  no  coating  system  case.  It  shows  the 


40 


smooth  variation  of  von  Mises  stress  distribution  along  the  depth  direction.  However,  for  the 
multilayered  FG  Ti/TiC/DLC  coating  system  case,  the  discrete  variation  of  von  Mises  stress 
profile  at  the  interfaces  exists.  Inside  the  DLC  coatings  and  FG  TiC  coatings,  higher  von  Mises 
stresses  were  generated  than  the  no  coating  case  because  of  the  higher  moduli  of  elasticity  for  the 
DLC  and  TiC  coating  materials.  If  more  distinct  layers  are  used  for  the  FG  multilayered  FEA 
model,  the  von  Mises  stress  profile  along  the  depth  direction  would  be  smoother,  that  is,  Ti  gs, 
Ti.75  and  Ti.65  are  used. 


DLC-2 


DLC-1 


Ti.25C.75 
Ti.30C.70 
Ti.50C.50 
Ti.70C.30 
Ti.90C.10 
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(b) 


Figure  4.9.  Effect  of  multilayered  FG  Ti/TiC/DLC  coating  system  on  von  Mises  stress 
distribution  for  the  indentation  depth,  8  =  50  nm,  and  coefficient  of  friction,  jj,  =  0.0. 


Figures  4.10(a)  through  4.10(d)  show  contour  plots  of  the  von  Mises  stress  distributions 
for  various  indenter  displacement  depths,  8  =  40  nm,  80  nm,  120  nm,  and  160  nm  at  the  end  of 
loading  stage.  The  indenter  displacement,  8  =  40  nm,  generates  the  maximum  Hertzian  pressure, 
Ph  =  4.5  GPa,  at  the  contact  for  no  coated  elastic  440C  stainless  steel  substrate.  The  indenter 
displacement,  8  =  160  nm,  generates  the  maximum  Hertzian  pressure,  Ph  =  9.0  GPa,  at  the 
contact  for  no  coated  elastic  440C  stainless  substrate.  Figure  4.10  (a)  depicts  the  contour  plot  of 
von  Mises  stress  distribution  inside  the  FG  multilayered  coating  system  generated  at  the  indenter 
displacement  depth,  8  =  40  nm.  The  maximum  von  Mises  stress,  3.64  GPa,  occurs  inside  the  FG 
TiC/DLC  coating  system.  Figure  4.10  (b)  displays  the  contour  plot  of  von  Mises  stress 
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distribution  inside  the  FG  multilayered  coating  system  generated  at  the  indenter  displacement 
depth,  8  =  80  nm.  The  increase  of  the  indenter  displacement  moves  the  maximum  location  of 
von  Mises  stress  further  below  the  contact  surface.  The  wider  contact  width  is  generated  due  to 
the  increase  of  the  indenter  displacement.  The  maximum  von  Mises  stress,  4.89  GPa,  occurs 
inside  the  FG  TiC/DLC  coating  system.  Figure  4.10  (c)  demonstrates  the  contour  plot  of  von 
Mises  stress  distribution  inside  FG  multilayered  coating  system  generated  at  the  indenter 
displacement  depth,  8  =  120  nm.  The  maximum  von  Mises  stress,  5.93  GPa,  occurs  inside  the 
FG  TiC/DLC  coating  system.  Figure  4.10  (d)  illustrates  the  contour  plot  of  von  Mises  stress 
distribution  inside  FG  multilayered  coating  system  generated  at  the  indenter  displacement,  8  = 

1 60  nm.  Due  to  the  magnitude  of  the  yield  strength  of  the  440C  stainless  steel  substrate,  plastic 
deformation  of  the  substrate  occurs.  The  maximum  von  Mises  stress,  6.57  GPa,  occurs  inside  the 
FG  TiC/DLC  coating  system.  Although  the  indenter  displacement  increases  from  8  =  40  nm  to  8 
=  160  nm,  the  location  of  maximum  von  Mises  stress  generated  at  the  contact  remains  inside 
DLC  coatings.  Also,  as  the  indenter  displacement  increases,  wider  contact  width  and  higher  von 
Mises  stress  occur. 
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Figure  4.10.  Contour  plots  of  von  Mises  stresses  distribution  generated  inside  the  multilayered  FG  Ti/TiC/DLC  coating 
system  from  the  FEA  model  at  various  indentation  depths;  (a)  5  =  40  nm,  (b)  8  =  80  nm,  (c)  8  =  120  nm,  and  (d)  8  =  160  nm. 
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Figures  4.1 1(a)  and  4.1 1(b)  show  the  von  Mises  stress  distribution  for  the  contact 

condition,  6  =  80  nm  and  p  =  0.0.  Figure  4.1 1(a)  depicts  the  von  Mises  stress  distribution 

generated  below  the  center  of  the  contact  along  the  depth  direction  on  the  plane  of  symmetry. 

Width-axis  corresponds  to  the  later  surface  distance.  Depth  =  0  representes  the  surface  position. 

This  indicates  that  the  maximum  von  Mises  stress  occurs  below  the  center  of  the  contact  inside 

DLC  multilayered  coating.  Due  to  the  distinct  material  properties  for  each  coating  layer,  the  von 

Mises  stress  distribution  displays  sharp  stress  changes  at  the  interface.  Figure  4. 1 1  (b)  illustrates 

the  contour  plot  of  von  Mises  stress  distribution  generated  inside  the  multilayered  FG 

Ti/TiC/DLC  coating  system. 

Max.  Von  Mises  Stress  occurs 

inside  Multilayered  Coating  x  10-6 


Figure  4.11.  (a)  Von  Mises  stress  distribution  on  the  plane  of  symmetry  generated  inside 
the  multilayed  FG  Ti/TiC/DLC  coating  system  for  the  indentation  depth,  8  =  80  nm  and 
coefficient  of  friction,  p  =  0.0  and  (b)  contour  plot  of  Von  Mises  stress  distribution. 

Figures  4.12(a)  through  4.12(d)  show  the  effect  of  applied  load  on  von  Mises  stress 
profiles  at  various  indentation  depths  for  the  contact  condition,  p  =  0.0.  Figure  4.12  (a)  depicts 
the  von  Mises  stress  profiles  along  the  depth  direction  at  the  center  of  the  contact  for  indentation 
depth,  8  =  40  nm,  for  both  FG  coated  system  and  440C  stainless  steel  substrate.  Without  any 
coating  on  the  440C  stainless  steel  substrate,  the  substrate  has  a  smooth  transition  of  von  Mises 
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stress  profile  along  the  depth  direction  with  maximum  2.37  GPa  located  below  the  surface. 
However,  the  mutilayered  FG  Ti/TiC/DLC  coating  system  has  a  drastic  transition  of  von  Mises 
stress  profile  along  the  depth  direction.  Especially  maximum  von  Mises  stress,  3.64  GPa,  occurs 
at  the  TiC/DLC  interface.  Also,  the  minimum  von  Mises  stress  occurs  at  the  bond  layer  with 
1 .46  GPa.  Figure  4.12  (b)  displays  the  von  Mises  stress  profiles  along  the  depth  direction  at  the 
center  of  the  contact  for  indentation  depth,  6  =  80  nm.  Without  any  coating  on  the  440C  stainless 
steel  substrate,  the  substrate  has  a  smooth  transition  of  von  Mises  stress  profile  along  the  depth 
direction  with  maximum  3.41  GPa  below  the  contact  surface.  Due  to  the  higher  moduli  of 
elasticity  for  TiC/DLC  materials  than  440C  stainless  steel,  the  von  Mises  stress  generated  inside 
the  FG  TiC/DLC  coating  system  has  much  higher  values  than  no  coating  case.  The  FG 
Ti/TiC/DLC  coating  system  has  a  drastic  transition  of  von  Mises  stress  along  the  depth  direction. 
Especially  maximum  von  Mises  stress,  4.89  GPa,  occurs  at  the  TiC/DLC  interface.  Also,  the 
minimum  von  Mises  stress  occurs  at  the  bond  layer  with  1 .46  GPa.  The  stress  generated  at  the 
bond  layer  interface  does  not  exceed  the  yield  strength  of  the  a  -  Ti  bond  layer  because  of  the 
elastic-perfectly  plastic  stress  strain  relationship  assumption.  Figure  4.12(c)  shows  the  von 
Mises  stress  profiles  along  the  depth  direction  at  the  center  of  the  contact  for  indentation  depth,  8 
=  120  nm.  Without  the  FG  Ti/TiC/DLC  multilayered  coating  system,  the  maximum  von  Mises 
stress,  4.00  GPa,  occurs  in  the  substrate  and  causes  the  plastic  deformation  inside  the  substrate. 
The  FG  TiC/DLC  coating  system  has  maximum  von  Mises  stress,  5.93  GPa,  at  the  DLC/TiC 
interface.  At  the  a-Ti  bond  layer,  the  von  Mises  stress  was  limited  to  1 .46  GPa  because  of  the 
material  properties  of  a-Ti  layer.  Figure  4.12(d)  depicts  the  von  Mises  stress  profiles  along  the 
depth  direction  at  the  center  of  the  contact  for  indentation  depth,  8  =  160  nm.  Without  the  FG 
Ti/TiC/DLC  multilayered  coating,  the  maximum  von  Mises  stress,  4.00  GPa,  occurs  in  the 
substrate  and  causes  the  larger  plastic  deformation  region  inside  the  substrate.  Although  there  is 
the  FG  Ti/TiC/DLC  multilayered  coating  system,  von  Mises  stress,  4.00  GPa,  occurs  inside  the 
substrate  and  causes  the  plastic  deformation  region  inside  the  substrate.  The  FG  TiC/DLC 
coating  system  has  maximum  von  Mises  stress,  6.57  GPa,  at  the  DLC/TiC  interface.  At  the  a-Ti 
bond  layer,  the  von  Mises  stress  was  limited  to  1 .46  GPa  because  of  the  material  properties  of  a- 
Ti  layer.  This  indicates  the 
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Figure  4.12.  Effect  of  applied  load  on  von  Mises  stresses  profiles  along  the  depth  direction  at  the  center  of  the  contact  at 
various  indentation  depths  for  friction  coefficient,  p  =  0.0;  (a)  8  =  40  nm,  (b)  8  =  80  nm,  (c)  8  =  120  nm,  and  (d)  8  =  160  nm. 


effect  of  bond  layer  on  the  von  Mises  stress  distribution,  that  is,  there  are  significant  changes  of 
von  Mises  stress  magnitude  at  the  bond  layer.  Generally  stress  profile  shape  remains  the  same. 
However,  the  magnitude  of  von  Mises  stress  increases  as  the  indenter  displacement  increases. 

Figures  4.13(a)  through  4.13(d)  show  the  contour  plots  of  equivalent  plastic  strain 
distribution  generated  for  FG  multilayered  coating  system  from  the  FEA  model  at  various 
indentation  depths.  Figure  4.13  (a)  illustrates  contour  plot  of  the  equivalent  plastic  strain 
distribution  generated  inside  the  FG  multilayered  coating  system  at  the  end  of  loading  stage.  The 
figure  indicates  that  the  maximum  equivalent  plastic  strain  occurs  at  the  a-Ti  bond  layer  for  the 
indenter  displacement,  8  =  40  nm.  The  plastic  strain  was  initiated  at  the  a-Ti  bond  layer  below 
the  center  of  the  contact  with  the  maximum  magnitude  of  9.875  x  10'  .  Figure  4.13  (b)  depicts 
the  contour  plot  of  equivalent  plastic  strain  distribution  generated  inside  the  FG  multilayered 
coating  system  for  the  indenter  displacement,  8  =  80  nm.  As  the  indenter  displacement  increases, 
the  magnitude  of  the  equivalent  plastic  strain  is  increased  to  6.133  x  10'2  at  the  a-Ti  bond  layer. 
The  location  of  the  maximum  equivalent  plastic  strain  moves  from  the  center  of  the  contact  to 
the  edges  of  the  contact.  Figure  4.13  (c)  demonstrates  the  contour  plot  of  equivalent  plastic 
strain  distribution  generated  inside  the  FG  multilayered  coating  system  for  the  indenter 
displacement,  8  =  120  nm.  As  the  indenter  displacement  increases,  the  magnitude  of  the 
equivalent  plastic  strain  is  increased  to  1.61 1  x  10'1.  The  location  of  the  maximum  equivalent 
plastic  strain  moves  further  from  the  axis  of  symmetry  to  the  edges  of  the  contact.  Figure  4. 13 
(d)  displays  the  contour  plot  of  equivalent  plastic  strain  distribution  generated  at  the  a-Ti  bond 
layer  for  indenter  displacement,  8  =  160  nm.  As  the  indenter  displacement  increases,  the 
magnitude  of  the  equivalent  plastic  strain  is  increased  to  2.598  x  10'1.  The  location  of  the 
maximum  equivalent  plastic  strain  moves  further  from  the  center  of  the  contact  to  the  edges  of 
the  contact,  i.e.  the  half  contact  width. 

During  the  nanoindentation  process,  the  multilayered  FG  Ti/TiC/DLC  coating  system  still 
experiences  the  residual  stress  after  the  nano  indenter  tip  is  unloaded  from  the  FG  multilayered 
coating  system.  Figures  4.14(a)  through  4.14(d)  show  the  contour  plots  of  residual  von  Mises 
stress  distribution  generated  inside  the  multilayered  FG  Ti/TiC/DLC  coating  system  obtained 
from  the  FEA  model.  As  the  indenter  tip  displacement  increases,  the  multilayered  FG 
Ti/TiC/DLC  coating  system  shows  the  higher  residual  von  Mises  stresses  generated  inside  the 
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Figure  4.13.  Contour  plots  of  equivalent  plastic  strain  distribution  generated  on  the  plane  of  symmetry  inside  the  multilayered  FG 
Ti/TiC/DLC  coating  system  at  various  indentation  depths;  (a)  8  =  40  nm,  (b)  8  =  80  nm,  (c)  8  =  120  nm,  and  (d)  8  =  160  nm. 
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Figure  4.14.  Contour  plots  of  residual  von  Mises  stress  distribution  generated  inside  multilayered  FG  Ti/TiC/DLC  coating  system  after 
the  unloading  of  nano  indenter  tip  from  various  indentation  depths;  (a)  5  =  40  nm,  (b)  8  =  80  nm,  (c)  5  =  120  nm,  and  (d)  8  =  160  nm. 
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coating  system.  Figure  4.14  (a)  depicts  the  contour  plot  of  residual  von  Mises  stress  distribution 
after  the  indenter  displacement,  5  =  40  nm,  is  applied  and  removed.  The  maximum  residual  von 
Mises  stress,  8.178  x  108  Pa,  is  generated  below  the  center  of  the  contact.  The  location  of  the 
maximum  residual  stress  is  at  the  a-Ti  bond  layer.  Figure  4.14  (b)  demonstrates  the  contour  plot 
of  residual  von  Mises  stress  distribution  inside  the  Ti/TiC/DLC  coating  system  after  the  indenter 
displacement,  5  =  80  nm,  is  applied  and  removed.  The  maximum  residual  von  Mises  stress 
generated  below  the  contact  is  1 .460  x  109  Pa  around  the  a-Ti  bond  layer.  As  the  indentation 
depth  increases,  the  larger  region  below  the  contact  surface  is  influenced  by  the  residual  stress. 
Figure  4.14  (c)  illustrates  the  contour  plot  of  residual  von  Mises  stress  distribution  inside  the 
Ti/TiC/DLC  coating  system  after  the  indenter  displacement,  5  =  120  nm,  is  applied  and  removed. 
The  maximum  von  Mises  stress  generated  below  the  contact  is  1.748  x  109  Pa  around  the  a-Ti 
bond  layer.  Figure  4.14  (d)  displays  the  contour  plot  of  residual  von  Mises  stress  distribution 
inside  the  Ti/TiC/DLC  coating  system  after  the  indenter  displacement,  5  =  160  nm,  is  applied  and 
removed.  The  maximum  residual  von  Mises  stress  generated  below  the  contact  is  2.244  x  109  Pa 
around  the  a-Ti  bond  layer.  Larger  area  of  the  coating  system  has  the  residual  von  Mises  stress 
inside  the  Ti/TiC/DLC  coating  system  after  the  indenter  tip  is  load  and  removed  from  the  contact 
for  the  indenter  displacement,  8  =  160  nm. 

Figure  4. 1 5  shows  the  effect  of  coating  system  design  on  von  Mises  stress  profiles  inside 
the  coating  systems  along  the  depth  direction  at  the  center  of  the  contact  on  the  plane  of 
symmetry  under  the  contact  condition,  8  =  200  nm  and  p  =  0.0.  The  dashdot  line  describes  the 
von  Mises  stress  profile  variation  from  the  substrate  surface  to  0.55  x  10'6  m  below  the  substrate 
surface  for  the  no  coated  substrate.  The  magnitude  of  the  von  Mises  stress  increases  upto  4  GPa 
as  the  depth  inside  the  substrate  increases  and  remains  constant  due  to  the  yield  strength  of  440C 
stainless  steel.  The  dashed  line  represents  the  von  Mises  stress  profile  variation  from  the  thin 
DLC  coating  surface  to  the  substrate.  The  maximum  von  Mises  stress,  5.67  GPa,  occurs  at  the 
DLC  coating  and  substrate  interface.  Also  the  maximum  von  Mises  stress  for  the  substrate  is 
limited  to  4.00  GPa.  The  dotted  line  reveals  the  stress  profile  variation  across  thickness  of  the 
thick  DLC  coating  and  substrate  system.  The  maximum  von  Mises  stress  of  10.1  GPa  occurs  at 
the  DLC  coating  and  substrate  interface.  Drastic  stress  profile  change  occurs  at  the  interface  due 
to  the  different  material  properties.  The  solid  line  shows  the  von  Mises  stress  profile  variation 


50 


generated  for  the  multilayered  FG  Ti/TiC/DLC  coating  system.  The  Ti/TiC/DLC  coating  system 
has  relatively  moderate  stress  variation  across  the  thickness  of  the  coating  system  compared  to 
the  thick  monolayered  DLC  coating  because  of  the  gradually  variation  of  material  properties 
across  the  thickness  of  the  coating  system. 
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Figure  4.15.  Effect  of  coating  system  design  on  von  Mises  stress  distribution  generated 
inside  the  coating  system  for  the  contact  condition,  §  =  200  nm  and  p  =  0.0. 

4.4  Summary  and  Conclusion 

Nanoindentation  technique  has  been  used  for  the  mechanical  characterization  of  coating 
materials.  The  finite  element  analysis  model  extends  the  usage  of  nanoindentation  technique  to 
the  calculation  of  material  property  such  as  yield  strength  from  the  experimental  nanoindentation 
data  such  as  hardness  and  modulus  of  elasticity.  Therefore,  nanoindentation  simulation  by  the 
finite  element  analysis  model  has  been  performed  to  determine  the  yield  strengths  for  the 
multilayered  FG  Ti/TiC/DLC  coating  materials  based  on  the  measured  nanoindentation  data. 
Using  the  measured  and  calculated  material  properties  for  the  coating  system,  the  FEA  model  has 
been  developed  to  investigate  the  stress  distribution  and  elastic-plastic  deformation  generated 
inside  the  multilayered  FG  Ti/TiC/DLC  coating  system  due  to  the  nanoindentation. 

The  superiority  of  the  multilayered  FG  Ti/TiC/DLC  coating  system  at  high  load 
applications  was  shown  and  discussed.  The  multilayered  FG  Ti/TiC/DLC  coating  system 
improves  the  system  integrity  under  the  heavy  loading  condition.  It  significant  reduces  the 
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magnitude  of  the  maximum  von  Mises  stress  generated  and  the  drastic  stress  profile  variation 
across  the  thickness  of  the  coating  system  than  monolayered  DLC  coating  system  due  to  the 
moduli  difference  of  the  coating  materials.  Also,  it  improves  the  wear  resistance  and  load 
support  mechanism  compared  to  both  no  coating  system  and  monolayered  DLC  coating  system. 
To  improve  the  adhesion  between  the  DLC/TiC  coating  system  and  substrate,  the  bond  layer  is 
used.  However,  due  to  the  material  properties  of  the  bond  layer  materials,  bond  layers  become 
the  prime  location  of  the  plastic  deformation  under  the  severe  high  loading  condition.  The 
residual  stress  should  be  considered  as  a  key  design  factor  when  the  multilayered  FG  coating 
system  is  developed. 
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